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Abstract 

The study continues the previous development [MATCH, 72 (2014) 39-73] of the per¬ 
turbative approach to relative stabilities of pi-electron systems of conjugated hydrocarbons 
modeled as sets of weakly-interacting initially-double (C=C) bonds. Distinct isomers of 
acyclic hydrocarbons (polyenes) are now under focus. The relevant total pi-electron en¬ 
ergies {£) are expressed in the form of power series containing members (£(k)) of even 
orders (k = 0,2,4,...) with respect to the averaged resonance parameter of initially-single 
(C-C) bonds. Terms to within the sixth order {k = 6) inclusive are shown to be of impor¬ 
tance for discrimination between similar isomers. In this connection, missing expressions 
for corrections £( 6 ) are originally derived. Conjugated paths of various lengths (i.e. linear 
chains consisting of C=C and C-C bonds alternately) are shown to be the most important 
(but not the only) fragments contributing to stabilization of any acyclic pi-electron system. 
Again, new types of fragments (substructures) are revealed (viz. the so-called composite 
conjugated paths) that contribute to destabilization of the system concerned. As a result, 
formation of the total energy of an acyclic pi-electron system is concluded to be governed 
by an interplay between stabilizing and destabilizing factors. Accordingly, the perturba¬ 
tive approach applied offers us an extension of the concept of conjugated paths. Particular 
isomers containing four, five and six C=C bonds are considered in a detail as examples. 


1 Introduction 

Qualitative intuition-based concepts and models play a crucial role in chemistry 
throughout its history. Accordingly, attempts to derive them deductively from more 


sophisticated quantum-chemical approaches contribute to our understanding of the 
nature of the given concept and/or model, as well as indicate directions for its possible 
extensions and improvements. 

The concepts of conjugated paths [1] and circuits [2-4] are successfully applied to 
evaluate relative stabilities of pi-electron systems of similar conjugated hydrocarbons 
[5], e.g. of various isomers of polyenes and of different Kekule valence structures 
of a certain benzenoid, respectively. Chains and cycles consisting of double (C=C) 
and single (C-C) bonds alternately are regarded here as the principal substructures 
determining stabilities of the structures concerned. Some limitations of these concepts 
also have been reported [5-8]. Difficulties in discriminating between stabilities of 
isomers of extended polyenes characterized by slightly different types of branching [5] 
are especially noteworthy in the context of the present study. 

In general, interpretation of stability of a certain pi-electron system depends on 
the model applied. Given that the latter coincides with the molecular graph [5, 9-11], 
the vertices and edges of which represent carbon atoms and carbon-carbon bonds, 
respectively, the relevant total energy is discussed in terms of properties of this graph 
(see e.g. [12,13]). Another alternative consists in modeling of a conjugated system as 
a set of weakly-interacting initially-double (C=C) bonds and thereby in employment 
of the perturbation theory to evaluate and to rationalize relative stabilities of isomers. 
Although the second option traditionally refers to acyclic conjugated hydrocarbons 
[14], adequacy of perturbative approaches to individual Kekule valence structures of 
benzenoids also is beyond any doubt [15]. Again, an evident analogy between the 
perturbational perspective to pi-electron systems and that underlying the concepts 
of conjugated paths and circuits gives us a hint about feasibility of a perturbation- 
theory- based derivation of these concepts followed by their extensions (if necessary). 

The above-formulated task, however, is not easily achievable. Difficulties in ap¬ 
plication to extended conjugated hydrocarbons of the most popular perturbational 
molecular orbital (PMO) theory [14] may be mentioned among the principal reasons 
here. Indeed, the original Dewar formula for the second order energy (£( 2 )) underlying 
this theory yields coinciding stabilities of numerous important pi-electron systems of 
distinct constitutions, including different Kekule valence structures of benzenoids and 



isomers of polyenes. To circumvent this difficulty, reference structures of artificial and 
involved nature are invoked instead of sets of C=C bonds, e.g. two allyle radicals for a 
Kekulc valence structure of benzene [14] . Such an option, however, makes the overall 
approach even more remote from the concepts of conjugated paths and circuits. 

An alternative to the standard PMO theory has been suggested recently [15,16], 
wherein corrections of higher orders (£( 3 ), £( 4 ), etc.) of the power series for total ener¬ 
gies (£) have been taken into consideration along with second order ones (£( 2 )) instead 
of passing to the above-discussed artificial reference structures. Thus, the classical 
model of conjugated hydrocarbons as sets of weakly-interacting C=C bonds has been 
preserved in this new approach. At the same time, the latter exhibited a much higher 
discriminative potential as compared to that of the standard PMO theory (which 
was shown to depend upon the number of corrections actually included). Besides, 
compact and chemically illustrative expressions for corrections £(3),£(4) [15-17] and 
£(5) [15] have been derived using an original matrix form of the perturbation the¬ 
ory (PT), namely the so-called non-commutative Rayleigh-Schrodinger perturbation 
theory (NCRSPT) [18-20]. Application of the latter approach to individual Kekule 
valence structures of benzenoid hydrocarbons has been carried out in the recent study 
[15]. Contributions both of linear (acyclic) and of cyclic conjugated fragments (sub¬ 
structures) were shown to be taken into consideration on the unified basis in the power 
series for total energies. This implies the approach employed to embrace perturbative 
analogues of both conjugated paths and circuits formally present in the structures 
concerned. The above-mentioned contributions, however, are not easily separable 
one from another. That is why no attempts were made to extend the qualitative 
concepts either of conjugated paths or of circuits in this study. 

The present article addresses a more simple case of acyclic conjugated hydro¬ 
carbons (polyenes) containing no conjugated circuits. Our aim now consists in a 
deductive revealing the principal types of conjugated fragments (substructures) con¬ 
tributing to the relevant total pi-electron energies and thereby in justifying and/or 
extending the concept of conjugated paths. At the same time, the extension being 
sought might be relevant also to numerous related concepts, such as conjugated paths 
invariants [ 21 ], the mean length of conjugated paths [ 22 ], conjugation paths used in 



studies of donor/acceptor functionalized pi-electron systems [23,24], etc. To achieve 
the above-formulated end, the same perturbative approach will be applied along with 
the experience of Ref. [15]. In respect of the latter, the following points deserve men¬ 
tioning: First, energy corrections of odd orders proved to take non-zero values for 
conjugated- circuits- containing systems only. Thus, we may now confine ourselves to 
members of the power series of even orders only. Second, energy increments £(o) and 
£( 2 ) were shown to be determined by total numbers of C=C and C-C bonds of the 
given structure, respectively. Thus, these corrections are expected to take coinciding 
values for isomers of the same hydrocarbon. In this connection, corrections at least 
of fourth and sixth orders seem to be required to discriminate between stabilities of 
these isomers. Thus, formulae for sixth order energies are originally derived in the 
Appendix of the present study. 

The paper starts with a brief overview of the principal expressions of the approach 
to be applied (Sect. 2). Thereupon, we turn to revealing the principal types of 
conjugated fragments that contribute to the energy corrections £( 4 ) and £( 6 ) (Sect. 
3). The final section is devoted to relative stabilities of specific polyenes. 

2 Expressions for energy corrections 

As already mentioned, acyclic conjugated hydrocarbons (polyenes) will be considered 
as sets of weakly-interacting initially-double (C=C) bonds. Moreover, the systems 
concerned belong to even alternant hydrocarbons (AHs) [9-11, 25,26]. The above- 
mentioned two points will be taken into consideration when constructing the relevant 
Hiickel type Hamiltonian matrix (H). 

Let the pi-electron system of a certain polyene to be initially represented by an 
2 N —dimensional basis set of 2 p z AOs of carbon atoms {y}, where N stands for the 
total number of C=C bonds. These AOs will be assumed to be characterized by 
uniform Coulomb parameters (a) as usual and the equality a = 0 will be accepted. 
As with the standard Hiickel model (see e.g. [9,25]), resonance parameters between 
AOs of chemically bound pairs of atoms only will be assumed to take non-zero values. 
Further, let the basis set {y} to be divided into two A-dimensional subsets {y*} and 



{x°} so that pairs of orbitals belonging to any chemical bond (C=C or C-C) find 
themselves in the different subsets. This implies the non-zero resonance parameters 
to take place in the off-diagonal (inter-subset) blocks of the Hamiltonian matrix (H). 
Accordingly, zero submatrices stand in the diagonal (intra-subset) positions of the 
matrix H as it is peculiar to AHs in general [9, 25,26]. Finally, let us enumerate 
the basis functions in such a way that orbitals belonging to the same C=C bond 
acquire the coupled numbers i and N + i. As a result, resonance parameters of these 
strong bonds take the diagonal positions in the intersubset blocks of the matrix H. 
Uniform values of these parameters (/?) also is among natural assumptions here. Let 
our (negative) energy unit to coincide with f3 in addition. The usual equality j3 — 1 
then immediately follows. Similarly, the averaged resonance parameter of weak (C-C) 
bonds will be denoted by 7 and supposed to be a first order term vs. the above- 
specified energy unit. 

In summary, Hamiltonian matrices of pi-electron systems of polyenes (H) take a 
common form that may be represented as a sum of zero (H(o)) and first order matrices 
(H(i)) including parameters of C=C and C-C bonds, respectively, viz. 


H 


H (o) + H (i) 


0 I 
I 0 


0 7 B 

7 B+ 0 ’ 


( 1 ) 


where I here and below stands for the unit matrix and the superscript + designates 
the transposed (Hcrmitian-conjugate) matrix. It deserves adding here that unit off- 
diagonal elements of the submatrix B ( B t] = 1 ,i / j ) correspond to C-C bonds, 
otherwise these take zero values. Meanwhile, the diagonal elements of the same sub- 
matrix (Bn) vanish because entire resonance parameters of C=C bonds are included 
into the zero order matrix H( 0 ). 

The Hamiltonian matrix of Eq.(l) coincides with that representing Kekule valence 
structures of benzenoids [15] because of similar constitutions of both systems. Ac¬ 
cordingly, the subsequent steps towards derivation of the relevant energy corrections 
also are similar. Thus, we start with passing to a new basis { 92 } consisting of bonding 
and antibonding orbitals of C=C bonds defined as normalized sums and differences 
of pairs of AOs \* and X°N+i an( l referred to below as bond orbitals. The transformed 
Hamiltonian matrix then meets the requirements of the NCRSPT (see the Appendix). 



As a result, general formulae for members of the power series for total energies (£) 
may be applied that have been derived earlier [15-17] using this PT [18-20]. As al¬ 
ready mentioned (Sect.l), we confine ourselves to terms £ p) of even orders (k = 0 , 2 , 
4,Let us turn now to individual members of the power series. 

The zero order energy (£( 0 )) coincides with 2N, whatever the specific structure of 
the given system. The subsequent second order member (£( 2 )) takes a rather simple 
form, viz. 

£( 2 ) = 4Tr(G ( 1 ) G+ ) ) > 0. (2) 

whereas the fourth order one (£( 4 )) consists of a sum of two components [16]: 

£<+> = 4Tr(G, 2 ) G+ ) ) > 0, £<“> = -4Tr(G ( 1 | G+ 1 ) G ( 1 | G+ ) ) < 0. (3) 

The notation TV here and below stands for a Trace of the whole matrix product 
within parentheses, and Gp) and G( 2 ) are the principal matrices of the NCRSPT 
[17-20] of the first and second orders, respectively, specified below. As is seen from 
Eqs. (2) and (3), Traces of positive-definite matrices [27] of the type AA + stand in 
these relations. Thus, sums of squares of elements of matrices Gp), G( 2 ) and Gp)G^ 
are contained there. This implies an a priori positive sign of the second order energy 
£( 2 ). Meanwhile, the components of the fourth order correction £( 4 ) are of opposite 
signs as indicated by additional superscripts (+) and (—). 

Let us now dwell on matrices Gp) and Gp) [17-20]. In the particular case of 
the NCRSPT employed in the present study (see the Appendix), these matrices are 
expressible as follows 

G,i) = -1r, G (2) = -1(SG (1) - G ( 1 ) Q) = 1(SR - RQ) = 1(SR + RS), (4) 

where matrices S, Q and R contain resonance parameters between the above-specified 
bond orbitals (BOs) { 93 }. Let bonding BOs (BBOs) and the antibonding ones (ABOs) 
to be correspondingly denoted by subscripts (+) and (—), e.g. <P{+)i and <f{-)i will 
stand for the BBO of the Ith C=C bond and for the ABO of the Lth one, respectively. 
Individual elements of matrices S, Q and R may be then explicitly expressed as follows 

£ij T(+)* | H | (^p_|_p >, Qi m <C ^(— )Z | H | )m >, Ril <C | H | >, 

( 5 ) 



where the BOs concerned are shown inside the bra- and ket-vectors. At the same 
time, the new matrices S, Q and R are related to the principal submatrices (7B and 
7 B + ) of our initial Hamiltonian of Eq.(l), viz. 

S = -Q=|(B + B + ), R=|(B+-B). ( 6 ) 

It is seen that matrices S(Q) and R are proportional to the symmetric (Hermitian) 
and skew-symmetric (skew-Hermitian) parts of the matrix B, respectively. On this 
basis, G(i) and G( 2 ) of Eq.(4) may be easily shown to be skew-symmetric (skew- 
Hermitian) matrices [28]. After an additional invoking the above-mentioned equality 
Bn = 0 for any i, we then obtain that 

'a Qii Rii G( 2 )m 0, (7) 

i.e. matrices embraced by Eq.(7) contain zero diagonal elements. A formal coincidence 
between matrices S and — Q also is seen from Eq.( 6 ). Just this circumstance allows 
us to eliminate the matric Q as shown in the last relation of Eq.(4). It deserves 
adding finally that the matrix product G(i)G^ determining the negative component 
of the fourth order energy (£^) is a symmetric (Hermitian) matrix. Consequently, 
diagonal elements (G(i)G^)jj take non-zero values and prove to be responsible for a 
large part of this energy component [16]. 

The energy correction of the sixth order (£( 6 )) is derived in the Appendix. Four 
components reveal themselves in this correction, viz. 

47 = 4Tr(GJ 3) Gg) > 0, (8) 

47 = 8 rr(G< 1 ,G+G( 1 ) G$ ) G ( 1 ) G+ ) ) > 0 , (9) 

47 = -32Tr[(G (1) G+ 2) )(G (1) G+ 2) ) + ] = -327V[(G+ ) G (2) )(G+ ) G (2) ) + ] < 0 , ( 10 ) 
4 ”) = 82>(G( 1 )G+ ) G (1| G+ ) ) = 8 rr[(G ( 1 | G+ ) )(G ( 2 | G+ ) ) + ], (11) 

where 

G(3) = - 1 (SG (2) - G (2) Q) = - 1 [(S) 2 R + 2SRS + R(S) 2 ). (12) 

The superscript o is used here to distinguish the above-introduced matrix G° 3 ^ from 
the standard third order matrix of the NCRSPT G^ defined by Eqs. (A5) and (A 6 ). 



The components of Eqs.(8) and (9) resemble of Eq.(3) in respect of both 
an a priori positive sign and skew-symmetric (skew-Hermitian) nature of underlying 
matrices G(C and GpjGiCGp), respectively. These components are correspondingly 
designated by additional subscripts 1 and 2. Accordingly, the only a priori negative 
component is shown in Eq.(10). The latter is alternatively expressible in terms of 
matrix products either G(pG^ or G)hG( 2 ). Meanwhile, the sign of the last com¬ 
ponent of the sixth order energy of Eq.(ll) cannot be established a priori and the 
superscript ( u ) (undefined) is used. 

Let us dwell now on interpretation of elements of the principal matrices deter¬ 
mining our energy increments of Eqs.(2) and (3) and (8)-(ll). Let us start with the 
simplest matrices Gp), G( 2 ) and G (A. As is seen from Eqs.(4) and (5), the element 
Gpyz connects the BBO (p(+)i and the ABO (p(-)i- Moreover, it is proportional to 
the relevant resonance parameter (Ru) and inversely proportional to the energy gap 
between BBOs and ABOs (equal to 2). Consequently, this element represents the 
direct (through-space) interaction between BOs (p(+p and <P(-)i- Besides, direct in¬ 
trabond interactions Gpy* vanish (see Eq.(7)). Again, the one-to-one correspondence 
between non-zero elements of the matrix B and C-C bonds along with Eq. (6) allows 
us to expect non-zero direct interactions (Gpyz ^ 0) to refer to BOs (<P(+)i and 
belonging to first-neighbouring C=C bonds only, the latter coinciding with those 
connected by a C-C bond. Further, the second order elements G^)u are accordingly 
interpretable as indirect (through-bond) interactions of the same BOs. Indeed, from 
Eq.(4) we obtain 


G(2)il , ^ ^ RimQml\i 


(13) 


(+)j (-) m 

where sums over (+)j and over (—)m correspondingly embrace all BBOs and all 
ABOs of the given system. It is seen that both BBOs (</?(+p) and ABOs (</?(_) m ) of 
other bonds play the role of mediators here [Note that j ^ i and m ^ l because of 
Eq.(7)]. Moreover, the orbitals <P(+)j and < P(-) m should overlap directly both with 
and with to be efficient mediators. That is why non-zero indirect interac¬ 
tions correspond to pairs of second-neighbouring C=C bonds possessing a common 
first neighbour. Analogously, the third order elements G° 3 y U represent the indirect 



interactions of the same BOs by means of two mediators. Pairs of BOs (ipt+u, (pi)) 
(</?i,</? 2 ) and {(p 2 ,(p(~)i) should overlap directly in this case, where ip\ and p >2 stand 
for mediating orbitals. 

Elements of matrix products determining the energy components £^,£^ 2 > £( 6 ^ 
and also may be interpreted as indirect interactions of BOs. For example, the 
element (G^G^hGp))^ represents a certain specific indirect interaction between BOs 
and wherein the mediating orbitals necessarily coincide with an ABO <P(-)\ 
and a BBO </?(+) 2 , respectively, whilst the interaction itself consists of three successive 
direct interactions G(i)u, Gm 12 and G(i) 2 ; referring to pairs of BOs (</?(+)*.,, <P(-)i) t 
(</?(_)!, </?( + ) 2 ) and (</?(+) 2 , Non-zero values of these direct components evidently 

are required to ensure a non-vanishing third order element (GpjG^Gp))^. Similarly, 
an element (G(i)G^)y involves a direct and an indirect interaction. Besides, pairs of 
bonding BOs play the role of interacting orbitals for elements both (G(i)G^Ojj and 
(G(i)Ge.g. the element (G(i)GiGjj represents the indirect interaction between 
BBOs <^( + )j and <P(+)j via ABOs of the first-neighbouring C=C bonds. It is evident 
that the mediating ABO should overlap with both <P(+)i and p>(+)j in this case too. 
Accordingly, the diagonal element (G(pG^)jj may be interpreted as the indirect 
self-interaction of the BBO <p(+)i- 

In summary, the above analysis yields the following rule: Any matrix element of 
the kth order connecting two BOs (p s and <p t and determining an energy component 
takes a non-zero value, if there is at least a single non-zero product of resonance 
parameters, i.e. 

< ip s I H I ip 1 >< ipi I H I ip 2 > ••• < (fk-2 | H I >< (p k -1 I H I (ft > 7 ^ 0, (14) 

where ip 1 , tp 2 , ...pk -1 stand for mediating orbitals. Given that the condition of Eq.(14) 
is met, we will say that in the given system there is a pathway of the (k-l)th order 
between BOs <p s and <pt- In the case of diagonal elements, we will accordingly have 
to deal with self-returning pathways. Besides, steps inside the same C=C bond are 
not allowed in these pathways because of Eq.(7). It also deserves emphasizing that 
the term a pathway (over BOs) is used here and below to make a distinction from 
conjugated paths defined in terms of chemical bonds. 



After returning to the power series for total energies of Eqs.(2),(3) and (8)-(ll), 
we may finally conclude that the higher is the order parameter (k), the more extended 
fragment of the whole system generally is embraced by the given correction (£(&)). 
In this respect, the present series resembles the graph-theoretic cluster expansion for 
total energy [29], as well as the expansion in terms of moments [30]. 

3 Conjugated fragments contributing to total en¬ 
ergies of polyenes 

As discussed already (Sect. 2) separate increments to total energies are determined 
by matrices Gp), G( 2 ), G^, G (1 )G^, etc. Thus, we will look for relations between 
elements of these matrices, on the one hand, and conjugated fragments present in 
the given system, on the other hand. The above-enumerated matrices are collected 
below into three groups that are analyzed separately. 

3.1 Relations between elements of matrices and G^\ and 
the simplest conjugated paths 

Let us start with elements of the first order matrix Gm defined by Eq.(4). An 
element G^u (as well as G^u) takes a non-zero value, if the Ith C=C bond and 
the Lth one (the underlying orbitals <P(+)i and belong to) are connected by a 
C-C bond (Sect. 2). This implies non-zero elements G^n(G^)u) to correspond to 
butadiene-like fragments and thereby to individual simplest conjugated paths (CPs) 
embracing two neighbouring C=C bonds and abbreviated below as CP(2)s. Moreover, 
the above-specified significant elements are local in their nature and, consequently, 
take uniform values for all CP(2)s. Let us also recall that the matrix Gm gives birth 
to the positive second order energy of Eq.(2). This implies all CP(2)s of the given 
polyene to contribute uniform stabilizing increments to the energy £( 2 ), the latter 
then being proportional to the total number of these paths. 

To exemplify the above simple rule, let us consider a linear polyene containing N 
C=C bonds and its cross-conjugated counterpart (dendralene) [31]. Carbon atoms 
and thereby the relevant 2p, AOs of these extended systems are assumed to be enu- 




I II III IV 


Figure 1: Isomers of octatetraene (I-IV) containing four C=C bonds (N=4). Num¬ 
berings of 2 p z AOs of carbon atoms also are shown, where AOs under numbers 1,2,3,4 
and 5,6,7,8 belong to subsets {y*} and {y°}, respectively. Any initially-double bond 
Ck—Cn+k is supposed to acquire the number K, where K=l,2,3,4. 

merated as follows 


_ p _ n (~a _ /-'i _ 

1 — ^N+l — ^2 — ^JV+2 — ^3 — ^AT+3--- — ~ '-'27V 7 

c, = c N+1 - c 2 (= Cat +2 ) - c N+3 (= c 3 ) - c 4 (= Cat+4) - 


where the 0=0 bonds placed outside the principal chain of carbon atoms of dendra- 
lene are shown within parentheses. Numbers 1,2...N and N+l, N+2...2N refer here 
to subsets {x*} and {y°}, respectively [see also Fig.l, where isomers of octatetraene 
I and IV serve as examples of the systems concerned for N=4], 

The principal first order matrices (viz. B and G(i)) of linear polyenes were shown 
to take a common form valid for any N [15] [these unified representations have be en 
denoted by B(IV) and G(i)(IV)]. After invoking Eqs.(l),(4) and (6), analogous matri¬ 
ces are easily constructable also for dendralenes. Let the latter to acquire additional 
superscripts /. For comparison, matrices G(i)(IV) and G 'm(N) are as follows 



0 

1 

0 

0 ... 


0 

1 

0 

0 ... 

G(,)(V)=-5 

-1 

0 

1 

0 ... 

. G ' wW—i 

-1 

0 

-1 0 ... 

0 

-1 

0 

1 ... 

0 

1 

0 

1 ... 


0 

0 

-1 

0 ... 


0 

0 

-1 0 ... 


where a standard factor (— 7 / 4 ) is introduced in front of matrices concerned for con¬ 
venience. It is seen that two non-zero elements (Gp)^ and Gp);*) correspond to any 
C-C bond in these matrices and thereby to any CP (2) of our polyenes, and these 
elements take uniform absolute values in addition. Moreover, matrices Gp)(IV) and 
Gp)(V) of Eq.(15) resemble one another except for signs of some elements. As a re¬ 
sult, the second order energies £^(N) and £'^(N) also are uniform and proportional 






to N — 1, viz. these coincide with 7 2 (N — l)/2. Such a result causes no surprise, as 
both linear and cross-conjugated polyenes contain the same numbers of C-C bonds 
and thereby of CP(2)s for the same N value. 

Let us turn now to elements of the second order matrix G( 2 ) defined by Eqs.(4) 
and/or (13). To ensure a non-zero value of the element G^)u (and of G ( 2 )u), the un¬ 
derlying BOs (p(+)i and <P(-)i should belong to second-neighboring C=C bonds (Sect. 
2), i.e. the Ith C=C bond and the Lth one should possess a common Erst neighbour 
coinciding with, say, the Mth C=C bond. The mutual arrangement of the three in¬ 
volved C=C bonds also playes an important role here: Given that the whole fragment 
1-M-L is of a linear constitution (see e.g. the isomer 1 of Fig.l), the mediating effect 
of the BBO <P(+) m and that of the ABO are added together and, consequently, 

a non-zero element G( 2 w(G( 2 m) arises. Meanwhile, the analogous increments cancel 
out one another for the cross-conjugated arrangement of C=C bonds I, M and J (e.g. 
in the isomer IV of Fig.l), and the relevant element G( 2 )u{G( 2 )u) vanishes. In sum¬ 
mary, two non-zero elements {G^)u and G( 2 )zi) correspond to BOs of the terminal 
C=C bonds (I and L) of any linear hexatriene-like fragment I-M-L and thereby to 
any conjugated path embracing three C=C bonds and further abbreviated as CP(3) 
[The remaining elements G( 2 ); m and G( 2 ) m ; vanish because of zero intrabond reso¬ 
nance parameters (see Eq.(7)]. This implies the total number of non-zero elements of 
a certain matrix G( 2 ) to coincide with the two-fold number of CP(3)s in the system 
under consideration. Uniform absolute values of the above-specified elements also 
easily follow from the definition of Eq.(4). As is seen from Eq.(3), the matrix G( 2 ) 
determines the positive (stabilizing) component (£^) of the fourth order energy that 
is an additive function with respect to squares of separate elements G( 2 )u(G( 2 )u) in 
addition. Consequently, the component consists of a sum of transferable incre¬ 
ments of individual CP(3)s and thereby it is expected to be proportional to the total 
number of these paths. 

The above-specified linear and cross-conjugated polyenes may be taken here again 
as examples. For the linear isomer, Eqs. (3), (4) and (15) yield the following common 



formulae 


g ( 2)(ao 


0 

0 

2 

0 

0 .. 

0 

0 

0 

2 

0 .. 

-2 

0 

0 

0 

2 .. 

0 : 

-2 

0 : 

0 : 

0 : 



87 4 (JV - 2) 
64 


(16) 


where 7 4 /64 is used here and below as a ’’subsidiary” unit of the fourth order en- 
ergy[15,16]. Elements of the matrix G( 2 )(iV) are chosen to coincide with 2 by choice 
of the front factor 7 2 /16 instead of 7 2 /8 [15,16] in order to reflect participation of 
mediating orbitals in pairs (e.g. (p(+) m and <p(-) m ) more conveniently. The expression 
for G (2) (N) of Eq.(16) illustrates the above-concluded one-to-one correspondence be¬ 
tween non-zero elements of the matrix G( 2 ) and individual CP(3)s. Proportionality 
between £^(N) and the total number of these paths (N — 2) also is seen. The 
fact that both G( 2 )(2) and £^ (2) vanish for butadiene (N=2) containing no CP(3)s 
causes no surprise here. By contrast, the alternating signs of elements when passing 
from one line of the matrix G^(A^) of Eq.(15) to another gives birth to a zero matrix 
GUN) for dendralenes in accordance with absence of CP(3)s in these hydrocar¬ 
bons. As a result, the stabilizing component of the fourth order energy S^'(N) also 
vanishes. 

Proportionality between the number of non-zero elements of the matrix G( 2 ) and 
that of CP(3)s deserves more illustration. To this end, let us consider the four isomers 
of octatetraene I-IV (Fig. 1). The matrix G( 2 )(/) and the energy increment £^ (I ) 
of the linear system 1 result directly from Eq.(16) for N=4, and £^(I) equals to 
167 4 /64. The remaining formulae under our interest are as follows 
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(17) 


Thus, total numbers of non-zero elements of matrices G( 2 ) coincide with two-fold 
numbers of CP(3)s for these systems too, i.e. with 4, 4, 2 and 0 for isomers 1-1V, 
respectively. Moreover, the related energy components {£[t) ) a l so are proportional 
to the same numbers. 






Therefore, a simple relation may be concluded between elements of matrices Gp) 
and G( 2 ), on the one hand, and the conjugated paths CP(2) and CP(3), on the other 
hand. Moreover, these CPs are the only conjugated fragments participating in the 
formation of elements concerned. Additivity of the consequent energetic increments 
8 ( 2 ) and with respect to transferable contributions of CP(2)s and CP(3)s, respec¬ 
tively, also is among conclusions here. 

Such a simple state of things, however, is no longer preserved when passing to 
terms of higher orders. To demonstrate this, we are about to consider elements of the 
matrix G" 3 ^ separately. 

3.2 Analysis of elements of the third order matrix G°^ 


Let us start with the above-discussed linear polyene containing N — 3 CP(4)s, where 
N> 3. The respective common third order matrix G° 3 ^(N) and the related sixth order 
energy increment (N) follow from Eqs.(8), (12) and (16), viz. 


G(3 ){N) = 



0 1 0 2 0 0 

-1 0 2 0 2 0 

0 -2 0 2 0 2 

-2 0 -2 0 2 0 

0 -2 0 -202 



4 7 6 [4(lV-3) + l] 
256 


(18) 


where — 7 3 /32 serves here and below as the standard factor for matrices G(A. Ac¬ 
cordingly, 7 6 /256 will be used as the ”subsidiary” sixth order energy unit. The first 
representatives of the series of matrices G° 3 ^(N) and of energy increments (N) 
corresponding to N=3 and 4 also deserve exhibiting, viz. 
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These particular cases evidently represent linear isomers of hexatriene and of oc- 
tatetraene I, respectively [Note that the energy increment £^1(3) follows directly 
from Eq.(18) after substituting N=3, but it is not the case for G° 3 -)(3). A more de¬ 
tailed discussion of this point may be found in Ref. [15]]. As is seen from Eq.(18), 






the energy increment £^(iV) contains a dependence upon the number of CP(4)s 
of the given chain (TV — 3) in accordance with the expectation. The total number 
of non-zero elements of the matrix G" 3 ^(TV), however, exceeds the two-fold number 
of CP(4)s considerably. Moreover, significant elements correspond not only to BOs 
of third-neighbouring C=C bonds (as it may be expected on the basis of the above 
experience), but also to orbitals of their first-neighbouring pairs. For example, the 
matrix G° 3 -)(4)[G^(/)] contains eight non-zero elements referring to BOs of C=C 
bonds under numbers (1,2), (1,4), (2,3) and (3,4) (Fig.l), and the elements concerned 
take non-uniform values in addition, e. g. 2 and 1 for pairs of BOs (+)1, (—)4 and 
(+)1, (—)2, respectively. It is evident that all the above-enumerated elements con¬ 
tribute to the stabilizing increment £^j(TV). Finally, neither the matrix G° 3 ^(3) itself 
nor the relevant energy correction £^j(3) vanish for the three-membered system of 
hexatriene (TV = 3) containing no CP(4)s [in contrast to the zero matrix G( 2 )( 2 ) 
and the vanishing correction 8^(2) of butadiene (TV = 2) discussed in the previous 
subsection]. 

To clarify the origin of these distinctive results, let us write down explicit expres¬ 
sions for elements of the matrix G"^(4)[G(' 3 ^(/)]. Using the first relation of Eq.(12), 
we obtain 

G(3)l 4 (-0 = —12 G(2)2 4 - G(2)13Q3 4 ), 

G ( 3 ) 23(-0 = - 2^ 21< ^ (2 ) 13 ~ G ( 2 ) 2 4 Q 4 3 ), 

G(3)i 2 (-0 = ^^i 2 ) 13 ^ 32 ’ G(3)34 (4) = --S 32 G(2)2 4 - (20) 

As is seen from the first formula of Eq.(20), a simple linear pathway from the BBO 
of the bond Ci=Cs (<^(+)i) to the ABO of the C^Cg (</?(-) 4 ) underlies the element 
G( 3 )i 4 (-0- Moreover, mediating effects of intervening orbitals (<P (+)2 and </?(_) 3 ) may be 
easily shown to be added together here. As a result, the absolute value of the element 
concerned coincides with 2. Thus, a simple relation may be concluded immediately 
between the element G" 3 - )14 (/) and the only CP(4) of the system I. A similar addition 
of contributing components takes place in the expression for the element G° 3 j 93 (J) too 
and the resulting value coincides with that of G° 3 ^ 14 (J). The underlying pathways, 
however, differ from linear ones in the latter case. Indeed, pathways over BOs con- 



taining self-returning segments correspond to both components of the expression for 
G( 3 ) 2,3 (I) > wherein orbitals of terminal bonds (<^(+)i and <P(-) 4 ) participate as media¬ 
tors. Meanwhile, orbitals of the fourth (C 4 =Cg) and first (Ci=C 5 ) bonds play no role 
in the formation of the first and second component of the element G^ 23 (/), respec¬ 
tively. It is evident that conjugated paths CP(3) embracing triplets of C=C bonds 
under numbers 1,2,3 and 2,3,4 may be correspondingly ascribed to the above-specified 
components. The same refers also to elements G° 3 ^ 19 (J) and G° 3 ^ 34 (J). 

We may expect, therefore, that matrices G° 3 ^ generally contain information not 
only about CP(4)s of the given polyene, but also about shorter conjugated paths. It 
is also likely that the presence of a standard CP (4) is not among necessary conditions 
for a non-zero matrix G" 3 j to represent a certain pi-clectron system. To support these 
anticipations, let us consider some polyenes of more involved constitutions. 

Let us start with the isomer of octatetraene III (Fig. 1) containing both conjugated 
and cross-conjugated fragments. As opposed to its linear counterpart I, the new 
isomer III contains no CP(4). Nevertheless, it is characterized by a non-zero matrix 
G" 3 j(//J) and a significant energy increment S^\(III), viz. 
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spite of the fact that the terminal C=C bonds (Ci=Cs and C 4 =Cs) are not joined 
by a conjugated path. This result causes little surprise if we recall the definition of 
the matrix G^ in terms of G( 2 ) shown in Eq.(12). Indeed, this definition indicates 
non-zero values of elements S 12 and G( 2)24 to be sufficient to ensure a significant 
element G " 3 ) 14 and this condition is met by the terminal orbitals and of 
our system III. Again, some similarity is beyond any doubt between constitutions of 
matrices G° 3 ^(I) of Eq.(19) and Gof Eq.(21). These two points allow the 
isomer III to be considered as a partially conjugated system. In this connection, a 
new concept of the semi-conjugated path may be introduced that contributes to the 
sixth order stabilization of the system along with the usual CP(4)s. In the present 
case, we will have to do with a semi-conjugated path embracing four C=C bonds 



and abbreviated below by SCP(4). The lower stabilizing effect of this new path as 
compared to the standard increment of the only CP (4) of the linear isomer I (see 
Eq.(19)) causes no surprise here. 

Another important distinction between the third order matrices along with 
the related energy increments and their counterparts of lower orders (Subsect. 
3.1) consists in much more involved dependences of the third (sixth) order char¬ 
acteristics upon the numbers of the standard conjugated paths when passing from 
linear to branched polyenes. This distinction may be traced back to the fact that the 
side subchains of branched systems offer new self-returning segments for pathways 
underlying separate elements G" 3 ^(G° 3 ^) and thereby the absolute values of these 
elements become excessively increased. For example, the branched octatetraene II is 
characterized by the following matrix Gand the energy increment £^j(/J) : 


G( 3 ) i J I) 



0 2 0 0 
-2011 
0-100 
0-100 



12y 6 

~ 256 ~ 


( 22 ) 


It is seen that absolute values of elements G° 3 ^ 12 (/J) and G" 3 ^ 21 (/J) coincide with 
G( 3 )u(^) and G° 3 ^ 23 (/) of Eq.(19) in spite of absence of CP(4)s in the branched 
system II. These increased elements are unambiguosly related to emergence of two 
self-returning segments in the pathways over BOs underlying elements G°^ 12 (II) and 
G( 3 ) 2 i (-f-0, namely of segments embracing the bonds C 4 =Cs and C 3 =C 7 . Neverthe¬ 
less, the total value of the stabilizing increment £^(//) is almost two times smaller 
for the branched isomer II as compared to the relevant value (207 6 /256) for its linear 
counterpart I. 

The decisive role of the side subchains in the formation of matrices G"^ (and 
thereby of energy components ) may be further illustrated by comparing these 
characteristics for isomers of decapentaene V, VI and IX (Fig. 2), where N=5. The 
total numbers of CP(4)s correspondingly equal to 2, 2 and 1 in these systems. Again, 
the isomers concerned may be regarded as consisting of the principal linear chain and 
of a side subchain, embracing the C=C bonds under numbers 1-4 and 5, respectively. 
The side subchain (C, 5 =Cio) then takes distinct positions with respect to the principal 
chain in the systems under comparison and thereby it offers different sets of self- 




IX X XI 


Figure 2: Isomers of decapentaene (V-XI) containing five C=C bonds (N=5). Num¬ 
berings of these bonds also are shown 


returning segments for pathways over BOs underlying particular elements of matrices 
G( 3 )(^ 0 > G? 3 )(V/) and G ° 3 ^(IX). As a result, distinctions may be anticipated both 
in the constitutions of the above-enumerated matrices and in values of the consequent 
energy components. 

To demonstrate this, let us start with matrices GjV. The matrix G^(I/)[G" 3 ^(5)] 
easily results from the general expression of Eq.(18) and contains eight elements equal 
to either 2 or -2 and four elements coinciding with either 1 or -1. Meanwhile, the 
remaining matrices and the consequent energy increments are as follows 
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It is seen that the matrix G° 3 ^(I/J) contains elements of higher absolute values as 
compared to G" 3 ^(/X) and £(g)j (VI) accordingly exceeds £^\ (IX). This result is 
in line with distinct numbers of CP(4)s in the isomers VI and IX (2 and 1) and 
thereby causes little surprise. It deserves adding, however, that the energy component 
(IX) does not coincide with £^\ (I) of Eq.(19) in spite of the same number of 
CP(4)s in the systems I and IX (equal to 1). Comparison of matrices G° 3 ^(V) and 
G ° 3 JVI) yields even more unexpected conclusions. Indeed, the matrix G ° 3 JVI) 



differs from G° 3 )(V) significantly in spite of the same number of CP(4)s present in 
both systems. Moreover, the sum of squares of elements of the former matrix exceeds 
that of the latter and, consequently, the branched isomer VI proves to be described 
by a higher energy increment £^\ as compared to its linear counterpart (V). This 
result implies a dependence of the sixth order energy upon the mutual arrangement 
of the two CP(4)s [Note that an increased overall stability of the branched system VI 
vs. the linear one (V) is not among the implications as discussed in Sect. 4], 

In summary, a relation is beyond any doubt between the third order matrix G° 3 ^ 
representing a certain polyene and the number of CP(4)s present there as it was 
the case with matrices Gq) and G( 2 ) determined by numbers of CP(2)s and CP(3)s, 
respectively (Subsect. 3.1). In contrast to the latter cases, however, the total number 
of CP(4)s is not the only factor determining the given matrix G^ and thereby the 
consequent energy increment : Other details of constitution of the given system 
also play their role here, e.g. presence of semi-conjugated paths (SCP(4)s) and a 
particular mutual arrangement of several CP(4)s (if any). 

Let us turn now to elements of matrix products of Eqs.(3) and (9)—(11). 

3.3 The role of composite conjugated paths in the formation 
of elements of matrix products 

Let us start with the product Gq)G^ determining the fourth order destabilizing 
energy component of Eq.(3). It is evident that an element (Gq)G(^)q does 
not vanish if in the given system there is an ABO such that Gqy; ^ 0 and 

G^y = G(i)_j 7 7 ^ 0. Because of the equality Gq)& = 0 for any i (see Eq.(7)), ABOs 
<P(-)i and <P(~)j are not able to play this role. Thus, the ABO <P(-)i necessarily belongs 
to a third (say Lth) C=C bond, where L^I and L^J. In other words, two simple 
mutually connected pathways are required here, namely a pathway from <P(+)i to 
c p(_)i and that from y?(-y and <P(+)j- Since a CP(2) corresponds to any element Gq^ 
(Subsect. 3.1), the above condition resolves itself into a requirement of two simple 
CP(2)s embracing a common (Lth) 0=0 bond. Given that this is the case, we will say 
that in the given system there is a composite conjugated path over three 0=0 bonds 
abbreviated below by COP(3). Besides, a COP(3) is automatically ensured under 



presence of a standard CP(3), but not vice versa (see Sect. 4). Further, the matrix 
product G(i)G(^ is characterized by non-zero diagonal elements (G(i)GtCjj that are 
interpretable as indirect self-interactions of respective BBOs tp(+)i via ABOs of the 
neighbouring C=C bonds (Sect. 2). In this connection, let us also define self-returning 
composite conjugated paths coinciding with squares of CP(2)s. As with the latter, 
these new paths also embrace pairs of neighbouring C=C bonds. Thus, let us use 
the abbreviation SRCCP(2). On the whole, the matrix G(i)G^ of a certain polyene 
may be then expected to contain information about both CCP(3)s and SRCCP(2)s. 

For illustration, let us consider matrices G(i)G^ representing the linear octate- 
traene (I), as well as its branched and cross-conjugated isomers II and IV (Fig. 1). 
The first two matrices under comparison differ one from another significantly, espe¬ 
cially in respect of diagonal elements, viz. 
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The reason for this distinction consists in the increased number of first neighbours of 
the second (C 2 =Cg) bond in the branched isomer II and thereby in the larger indirect 
self-interaction of the relevant BBO v 7 (+) 2 - In terms of conjugated paths we have to 
do here with an increased number of SRCCP(2)s referring to the 2nd C=C bond. 
If we recall here that isomers I and II both contain two CP(3)s, the above result 
becomes even more important in distinguishing between their stabilities (Sect. 4). So 
far as the matrix G(i)G^(/V) is concerned, it resembles G^GtJl) in respect of 
absolute values of all matrix elements [negative off-diagonal elements of G(pG^(/) 
become replaced by positive ones when passing to G(i)G^(/V)]. This implies the 
linear isomer I and the cross-conjugated one IV to contain the same sets of both linear 
and self-returning composite conjugated paths (i.e. of CCP(3)s and SRCCP(2)s) in 
spite of different numbers of the standard CP(3)s (2 and 0). 

Let us turn now to third order matrix products and start with G(i)GtG(i). 
The skew-symmetric (skew-Hcrmitian) nature of this matrix (Sect. 2) implies zero 
values for any diagonal element (GpjG^Gp))^. Nevertheless, an analogy still exists 



between products GpjG^Gp) and G(i)G^. In particular, the matrix G(i)G^G(i) 
may be similarly shown to represent products of three connected CP(2)s, the total 
number of the embraced C=C bonds generally coinciding with four. Accordingly, we 
may define composite conjugated paths over four C=C bonds (CCP(4)s). It should be 
mentioned, however, that the overall situation becomes somewhat more involved when 
passing from G(i)G^ to G^G^Gp), as it was the case when comparing matrices 
G( 2 ) and G° 3 ) (Sect. 3). In particular, pathways over BOs underlying elements 
(G(i)G( l 1) G(i))^ may possess self-returning segments (as it was the case with G° 3 ^). 
Consequently, the resulting CCPs actually embrace three or even two C=C bonds. 
It is also evident that zero diagonal elements of the product GpjG^Gp) reflect 
impossibility of completely self-returning paths in this case. Finally, matrix products 
G(i)G(^ and G( 2 )G(^ remain to be discussed. In contrast to previous cases, these 
products are neither symmetric (Hermitian) nor skew-symmetric (skew-Hcrmitian) 
matrices. It is evident that the underlying CCPs consist of products of a CP(2) and 
a CP(3), and of a CP(3) and a CP(2), respectively, and also generally embrace four 
C=C bonds. Self-returning segments are possible here too. 

It is seen, therefore, that elements of matrix products are determined by conju¬ 
gated paths of a non-standard (viz. composite) nature. Moreover, most of these new 
paths give birth to destabilizing energy components (e.g. and £^)- 

Before finishing this Section, the following remark deserves to be made: Energy in¬ 
crements originating from matrix products G(pG^, G(i)G^, etc. and these related 
to simple matrices G(i),G( 2 ) and G" 3 ^ hardly are independent. Quite the reverse, a 
certain interdependence may be foreseen between some of these increments, e.g. be¬ 
tween and The main reason for such an anticipation consists in the presence 
of the same matrix G( 2 ) in the definitions of underlying matrices G° 3 ^ and G(i)G^ 
(see Eq.(12)). [Besides, both matrices G° 3 ^ and G(i)G(L vanish, if G( 2 ) coincides 
with a zero matrix]. Thus, an increased stabilization of a certain system due to a 
larger increment may be expected to be accompanied by a growing destabiliza¬ 
tion (£( 6 J) and vice versa. Given that relative stabilities of different isomers are under 
interest (Sect. 4), the outcome of the comparison then depends on the overall balance 
between energy increments of opposite signs. 



4 Discussions of relative stabilities of specific iso¬ 
mers 


Let us start with comparison of linear polyenes to their cross-conjugated isomers 
(dendralenes) in respect of overall relative stabilities of their pi-electron systems. 
Zero order energies of these isomers are uniform for the same number of C=C bonds 
(N), viz. £( 0 ) (N) = SUN) = 2 N. The relevant second order increments also coincide 
one with another (see the implications of Eq.(15)). Finally, the stabilizing compo¬ 
nents of the fourth order energies are discussed in the Subsect. 3.1 [S^(N) is shown 
in Eq.(16) and takes a significant value because of N-2 three-membered conjugated 
paths (CP(3)s) present in the linear system, whereas S^'(N) vanishes due to ab¬ 
sence of these paths in dendralenes]. Thus, let us turn immediately to the remaining 
(destabilizing) components of the fourth order energies. 

Employment of the expression for the matrix Gp )(N) of Eq.(15) to construct 
the product (G(i)G^)(iV) shows this important representation of linear polyenes to 
contain elements 1,2,2...2,1 in its principal diagonale along with —1 in the second- 
neighbouring off-diagonal positions [The matrix (G(i)GmX-O of Eq.(24) serves as 
an example for N=4], Moreover, the analogous matrix product (G(i)G^)'(V) rep¬ 
resenting the dendralene series and originating from G'^JN) of Eq.(15) also closely 
resembles the above-discussed one, except for opposite (i.e. positive) signs of all off- 
diagonal elements. Accordingly, destabilizing components of fourth order energies are 
uniform for both systems under comparison. The latter result may be entirely traced 
back to an evident fact that products of CP(2)s and thereby numbers of composite 
conjugated paths CCP(3)s coincide one with another in both polyenes. The relevant 
total fourth order energies are then as follows 
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where S^ (V) is taken from Eq. (16). It is seen that the correction £( 4 ) (V) is a positive 
quantity for N^4 owing to predominance of its stabilizing component £^\n) over 
the destabilizing one Sj^ (N). This implies the sufficiently long linear polyenes to 
be additionally stabilized vs. the sum S^(N) + £( 2 )(N) due to the fourth order 
energy. By contrast, the analogous correction of dendralenes (S'^(N)) consists of the 



destabilizing component only and, consequently, it is a negative quantity [Particular 
cases of Eq.(25) referring to N=2 and N=3 also are of interest. In the case of butadiene 
(N=2), both relations yield the same correction £( 4 ) (2) equal to —2y 4 /64 in accordance 
with the expectation, which is a negative quantity in addition due to absence of 
CP(3)s. This correction consists of the destabilizing component only that originates 
from SRCCP(2)s. For the linear hexatriene (N=3), the two components (£,^(3) and 
£( 4 )\s)) cancel out one another and the total fourth order energy takes a zero value. 
The branched isomer of hexatriene, in turn, is characterized by a negative fourth 
order energy]. 

Thus, a higher relative stability of linear polyenes vs. dendralenes unambiguosly 
follows from our results and this conclusion coincides with those of other approaches 
[32-35]. As with the standard model of conjugated paths [5], the above analysis also 
indicates the presence of CP(3)s to be the origin of the increased stability of linear 
isomers. In contrast to the standard model, however, an additional destabilizing 
factor is now revealed to manifest itself in both systems under comparison that is 
interpretable as a contribution of composite conjugated paths. Only because of the 
above-established coincidence of absolute values of the underlying energy increments 
for linear and cross-conjugated systems, the destabilizing factor becomes irrelevant 
when comparing their relative stabilities. 

For other series of polyenes (e.g. the branched ones), a general analysis like that 
carried out above hardly is feasible. Thus, we will confine ourselves to comparisons 
of relative stabilities of individual representatives of different series. 

Let us start with the four isomers of octatetraene I-IV (Fig.l) containing the same 
number of C-C bonds and thereby of CP(2)s (equal to 3). Accordingly, the relevant 
second order energies also are uniform. Again, the total numbers of CP(3)s corre¬ 
spondingly coincide with 2, 2, 1 and 0 for systems I-IV. Since the linear isomer (I) 
and its cross-conjugated analogue (IV) are particular cases of the above-considered 
polyenes, Eq.(25) yields £( 4 )(/) and £( 4 ){IV) equal to 2y 4 /64 and —147 4 /64, respec¬ 
tively. 

Let us now dwell on the branched isomer II. The stabilizing component of the 
fourth order energy £^\ll) is shown in Eq.(17) and coincides with £^ (/) following 



from Eq.(16) for N=4 owing to similar non-zero elements of matrices G( 2 )(//) and 
G( 2 )(/) (Subsect. 3.1). This result is in line with the same number of CP(3)s in 
polyenes I and II. Meanwhile, the matrix G(i)G^(//) differs from both G(i)GV(/) 
and G(i)G^(/ld) significantly (see Eq.(24) and the discussion nearby). It is evident 
that the sum of squares of elements of the matrix G(i)GV(//) exceeds the relevant 
value for isomers I and IV. As a result, the destabilizing component S^(II) is of an 
increased absolute value vs. £^\l) and 8 ^ {IV). The overall result referring to the 
isomer II is then as follows 
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Hence, the total fourth order energies £^{I) and £^(II) differ one from another 
for isomers I and II in spite of the same numbers of CP(3)s. Moreover, the linear 
isomer is predicted to be more stable as compared to the branched one (II). Although 
this result is in line with predictions of the standard CP model (as well as with the 
relevant general graph-theoretical results [34,35]), the above analysis indicates another 
underlying reason. Indeed, the branched isomer II is now concluded to be less stable 
owing to a greater destabilizing effect of self-returning composite conjugated paths 
(SRCCP(2)s) defined in Sect. 3. Meanwhile, the greater stability of the isomer I (vs. 
the branched analogue II) is traced back to the presence of a CP (4) in the linear chain 
when the usual model of conjugated paths is applied. 

Finally, the semi-conjugated isomer (III) remains to be discussed. The relevant 
matrix G( 2 )(///) is shown in Eq.(17) and contains two non-zero elements in ac¬ 
cordance with a single CP (3) present in the given system. Meanwhile, the matrix 
G(i)G( l jj(///) closely resembles G^GJM) of Eq.(24) in respect of absolute values 
of non-zero elements. We then obtain 
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Thus, the above results indicate the following order of relative stabilities of isomers: 
I > II > III > IV. Completely similar conclusions follow also for analogous isomers 
of decapentaene V, VI, VII and VIII shown in Fig. 2. The relevant fourth order 
energies correspondingly coincide with 4y 4 /64, 0, —4y 4 /64 and —12y 4 /64.Besides, 



the above-concluded relative orders of stability for both I-IV and V-VIII are in line 
with the maximal it— energy of 1,1-divinyl isomers of polyenes among the branched 
ones established in Ref. [35]. 

It is seen, therefore, that fourth order energies are sufficient to distinguish between 
relative stabilities of representatives of different principal series of polyenes. In the 
case of distinct isomers characterized by more similar overall constitutions, however, 
the fourth order energies often are uniform and, consequently, sixth order corrections 
should be invoked. Let turn now to relevant examples. 

Let us start with comparison of the above-discussed isomers of decapentaene VI 
and IX containing the same number of CP(3)s equal to three. As it may be easily 
proven after constructing the relevant principal matrices, a zero fourth order energy 
is peculiar to both isomers under comparison and this result evidently causes little 
surprise. Again, distinct numbers of CP(4)s (namely 2 and 1 for systems VI and IX, 
respectively) allow us to expect the sixth order energies to be responsible for different 
stabilities of these isomers (Subsect. 3.2). Thus, let us now turn to corrections 
S{%){VI) and S {6) (IX). 

Matrices G° 3 JVI) and G° 3 ^(IX) along with the consequent energy increments 
^(6)i (VI) and (IX) are shown in Eq.(23) and discussed nearby. A higher value 
of £(6 )i (VI) vs. (IX) was in line with the relevant numbers of CP(4)s. Further, 
matrices G(i)G ( t 1) G(i)(VJ) and G^G^Gp^/X) determining the second stabilizing 
increments S^ 2 {VI) and S^ 2 (IX), respectively, prove to be similar in respect of 
absolute values of their non-zero elements. This implies the numbers of the relevant 
composite conjugated paths (consisting of three CP(2)s) to be uniform in the isomers 
VI and IX. As a result, energy increments £^ 2 (VI) and £^\(IX) also coincide one 
with another and equal to 407 6 /256. 

As with the above-discussed matrices G? 3 )(V7) and G ° 3 ^(IX) of Eq.(23), the ma¬ 
trix GpjGpMV/) also contains more non-zero elements as compared to G(i)G^(/X) 
[Apart from a single element equal to 2, these matrices involve eight and seven el¬ 
ements, respectively, that are equal to either 1 or —1. This result is in line with 
the above-foreseen parallelism between alterations in elements of matrices G° 3 ) and 
G(i)G( 2 ) (see the last paragraph of Sect. 3)]. Consequently, the absolute value of the 



sixth order destabilizing increment 8^ (VI) also exceeds that of 8^ (IX). Moreover, 
the same refers also to increments 8^ (VI) and 8^ (IX), viz. 



(28) 


As is seen after summing up the relevant contributions, the sixth order stabilization 
energy increases by 187 6 /256, whereas the absolute value of destabilization grows 
only by 87 6 /256 when passing from IX to VI. The total sixth order energies then 
coincide with zero and —107 6 /256 for isomers VI and IX, respectively, and indicate 
the former pi-electron system to be more stable than the latter in accordance with 
graph-theoretical conclusions of Ref. [35]. The present result may be traced back 
to the relevant numbers of CP(4)s. The relation between these numbers and relative 
stabilities of isomers VI and IX, however, is far from being of a straightforward nature 
as the above discussion shows. 

Let us return again to the semi-conjugated isomer of decapentaene (VII) and 
compare it to a similar one (X) (Fig. 2). As with the above-considered couple (VI and 
IX), the isomers VII and X also are characterized by the same numbers of CP(3)s and, 
consequently, by uniform fourth order energies equal to —47 4 /64. Different numbers 
of CP(4)s of these hydrocarbons also deserve mention here (these coincide with 1 and 
0 for isomers VII and X, respectively). Thus, let us turn to the sixth order energies 
8 (6) {VII) and £ m (X). 

Matrices G °^(VII) and G° 3 ^(V) exhibit a clear parallelism between absolute 
values of their non-zero elements and the respective numbers of CP(4)s as previ¬ 
ously. Accordingly, the total number of non-zero elements is higher in the matrix 
G(i)G ~V(yiI) as compared to G(i)G^(X). Meanwhile, matrices G( 1 )G^G( 1 )(V//) 
and G(!)G^G(i)(V) contain analogous non-zero elements. The total sixth order en¬ 
ergies 8(6){VII) and 8^(X) then correspondingly equal to 2^/256 and —87 6 /256. 
Thus, the isomer VII is predicted to be more stable as compared to X in analogy with 
the above-considered couple VI and IX. The decisive role of CP(4)s in the formation 
of this result also is beyond any doubt. 

Let us now compare the semi-conjugated system VIII to a closely related one (XI). 
These isomers also are characterized by coinciding numbers of CP(3)s and by uniform 



fourth order energies (equal to —127 4 /64). As opposed to previous examples, how¬ 
ever, both VIII and XI contain no CP(4)s. This implies that the standard model of 
conjugated paths is not able to distinguish between their stabilities. Again, the same 
isomers VIII and XI differ one from another in numbers of semi-conjugated paths 
(SCP(4)s) defined in Sect. 3. Indeed, the system VIII contains a single SCP(4) em¬ 
bracing the C=C bonds under numbers 2,3,4,5. Meanwhile, the remaining isomer XI 
involves two SCP(4)s made np of C=C bonds 1,2,3,4 and 2,3,4,5. In this connection, 
comparison of sixth order energies of polyenes VIII and IX is of particular interest. 

As with the above-considered couples of isomers, matrices G(i)G^G(i)(V7i7) 
and G(i)G^G(i)(XJ) contain uniform sets of non-zero elements and contribute co¬ 
inciding increments to the relevant sixth order energies. Again, matrices G °^(VIII) 
and G° 3 }(XI) correspondingly involve six and eight non-zero elements equal to either 
1 or —1 and thereby reflect different numbers of SCP(4)s present in these systems. 
Consequently, the energy increments £^\(V 111) and £^(AV) coincide with 6^/256 
and 87 6 /256, respectively. Thus, £^(V/) exceeds (VIII) by 27 6 /256 in accor¬ 
dance with the expectation. Analogously, matrices Gp)G t^iVIII) and G(i)G^(A"J) 
contain three and four unit elements, respectively, and yield the following energy in¬ 
crements 





32y 6 

~256”’ 


8$ (VI11) 



0. (29) 


Thus, the absolute value of the destabilizing increment is increased by 87 6 /256 when 
passing from VIII to XI and this alteration exceeds that of the stabilizing increment 
considerably. It is no surprise in this connection that the isomer VIII proves to 
be more stable as compared to its counterpart XI [The total sixth order energies 
correspondingly equal to 107 6 /256 and 47 6 /256]. Such a somewhat unexpected result 
may be entirely traced back to the higher destabilizing effect of composite conjugated 
paths (CCP(4)s) underlying matrices G(i)G^ in the isomer XI vs. VIII. 

The last example under our interest embraces four isomers XII, XIII, XIV and XV 
of Fig. 3, all of them containing six C=C bonds (N=6) and characterized by uniform 
fourth order energies —67 4 /64 in addition. The total numbers of CP(3)s also are 
uniform here and coincide with 3. Meanwhile, the relevant numbers of the standard 



1 



XIV XV 


Figure 3: Selected polyenes (XII-XV) containing six C=C bonds (N=6). Numberings 
of these bonds also are shown 


CP(4)s correspondingly equal to 0, 1, 1, and 2. This implies the isomers XIII and 
XIV to be of the same composition in terms of the standard conjugated paths. [It is 
no surprise that these isomers have been never discriminated as concluded in Ref. [5]]. 
The relevant numbers of SCP(4)s coincide with 3,2,1 and 1, respectively. 

Separate increments to the sixth order energies of the above-enumerated systems 
also follow the above-observed trends. Thus, the first stabilizing increments ) 
originating from matrices G° 3 ^ are as follows 


sftliXII) = ^-,E)VAXIII) = ^-,S}VAXIV) = —,C(AT) = 

(6)iV ) 256 (6)lV ’ 256 (6)lV ’ 256 (6)lV ’ 


227 6 „ (+ ) 


34/ -(+), 


32 / -(+), 


42y 6 

~256 


(30) 

and correlate with the total numbers of both CP(4)s and SCP(4)s. Furthermore, ma¬ 
trices G(i)G/G(!) are similar as previously except for the matrix G(i)G/G(i )(XIII) 
possessing a somewhat higher number of non-zero elements as compared to the re¬ 
maining ones. The relevant energy increments the take the form 


= £[&(XIV) = £<+> (XV) = £$(xnil) = ^ 


76 ) 2 '' 


76)2' 


(31) 


Finally, matrices G(i)G^ contain different numbers of unit elements along with a 
single element 2, namely seven, nine, eight and nine unit elements for isomers XII- 



XIV, respectively. These matrices give birth to following energy increments 




£$( XIII ) = £$( XV ) = - 

0, £j$(XIII)=£${XIV) = ^, 



967 s 

~256”’ 

(32) 


It is seen that absolute values of destabilizing increments also correlate with total 
numbers of both CP(4)s and SCP(4)s in the isomers XII-XIV. After summing up all 
increments concerned, we obtain 

%,(aw) = -% e m (xnt) = s^xiv) = -g£, e^(xv) = ^ 

(33) 

Thus, the relative stability grows with the increasing number of the standard CP(4)s 
in this case too. Moreover, the isomers XIII and XIV (both containing a single CP(4)) 
also are discriminated when applying the perturbative approach: The isomer XIII is 
predicted to be less stable as compared to XIV as it was the case with XI vs. VIII. 


5 Conclusions 

Analysis of power series for total energies of pi-electron systems of acyclic conjugated 
hydrocarbons (polyenes) supports the principal assumptions underlying the model of 
conjugated paths and thereby offers a justification of the latter. In this respect, the 
following points may be mentioned: 

(i) The standard conjugated paths (CPs) embracing two, three and four linearly- 
connected C=C bonds (CP(2)s, CP(3)s and CP(4)s) contribute significantly to terms 
of power series for total energies of the second, fourth and sixth orders, respectively; 

(ii) The afore-mentioned contributions always are of positive signs (in negative 
energy units) and thereby of stabilizing nature; 

(iii) Relative values of the contributions concerned depend upon total numbers 
of the respective standard CPs present in the given hydrocarbon. In particular, 
contributions of CP(2)s and CP(3)s to energy corrections of the second and fourth 
orders, respectively, are expressible as sums of transferable increments of individual 
CPs and thereby these are directly proportional to the numbers of the latter; 



(iv) The decisive energy correction of the second order (coinciding with the De- 
war energy of the PMO theory [14]) is determined exclusively by the number of the 
simplest conjugated paths (CP(2)s) embracing two C=C bonds connected by a C-C 
bond. 

Again, application of the perturbative approach to relative stabilities of pi-electron 
systems of polyenes undertaken in the above study contributes to an extension of the 
very concept and/or model of conjugated paths. This conclusion is based on the 
following properties of the power series for total energies: 

(i) Members of the series of the fourth (£( 4 )) and sixth orders (£(e)) contain both 
positive (stabilizing) and negative (destabilizing) components and these are in some 
relation one with another in addition. This implies that destabilizing factors also 
manifest themselves in polyenes that are able to play an equally decisive role in the 
formation of the final total energy; 

(ii) Negative (destabilizing) components of the energy corrections £( 4 ) and £( 6 ) 
are interpretable as contributions of conjugated paths of a non-standard (composite) 
nature defined as successive products of two or three connected standard CP(2)s 
and/or CP(3)s. Moreover, the self-returning composite conjugated paths (SRCCPs) 
prove to be especially important, wherein an even number of C=C bonds is involved 
and each of them is visited twice; 

(iii) The usual linear conjugated paths embracing four C=C bonds (CP(4)s) are 
not the only fragments (substructures) participating in the formation of positive (sta¬ 
bilizing) components of the sixth order energy. In particular, four-membered frag¬ 
ments containing both a linearly-conjugated segment and a cross-conjugated one [the 
so-called semi-conjugated paths (SCP(4)s] also contribute to the sixth order stabi¬ 
lization of polyenes; 

(iv) The energy correction of the sixth order generally is a non-additive quantity 
with respect to increments of individual participating fragments (substructures) in¬ 
cluding the standard CPs. Consequently, the total value of this correction depends 
upon the actual mutual arrangement of these fragments in the given hydrocarbon. 

Consideration of specific examples also corroborates the extended nature and, 
consequently, a higher discriminative potential of the perturbative approach applied 



vs. the usual CP model. In this respect, the most important conclusions are as 
follows: 

(i) Destabilizing increments of self-returning composite conjugated paths (SRC- 
CPs) prove to be generally responsible for lower relative stabilities of branched isomers 
of polyenes vs. their linear counterparts; 

(ii) isomers of extended polyenes containing different numbers of semi-conjugated 
paths (SCPs) usually are represented by distinct sixth order energies even if the 
relevent numbers of the standard CPs are uniform. 

Appendix 

A Derivation of expressions for the sixth order en¬ 
ergy corrections 

In its most general form, the power series for total energies of molecules and molecular 
systems has been originally derived in Ref. [17]. This study contains members of the 
power series up to fourth order (k = 4). The relevant fifth order terms may be 
found in Ref. [15]. A direct extension of the above-cited derivation to terms of higher 
orders (including the sixth order ones) is a rather cumbersome procedure. In this 
connection, we will confine ourselves here to a less general Hamiltonian matrix vs. 
that of Refs. [17,18] as described below. Nevertheless, the overall methodology to 
be invoked closely resembles the original one [17]. The main points of the latter are 
as follows: First, the interrelation [36] is employed between the total energy being 
sought (£), the Hamiltonian matrix of the system(s) concerned (H) and the relevant 
representation of the one-electron density matrix (the charge- bond order (CBO) 
matrix) P, viz. 

S = TV (PH). (Al) 

Second, the matrix P is derived directly [18] on the basis of solution of the so-called 
commutation equation [36]. For Hamiltonian matrices (H) consisting of zero and 
first order members (H( 0 ) and H^), respectively), the above-mentioned solution may 
be carried out perturbatively. As a result, both the CBO matrix P and the total 
energy S are expressible as sums of corrections P(^ and of increasing orders (k). 



Moreover, each energy correction Sn.) is additionally representable as a sum of two 
components, viz. 


%l = 4“) + £ W • 4t = rr(P w H (0) ), £% = TrfP^DH,,,). (A 2 ) 


The most general Hamiltonian matrix (H) underlying the original derivation of 
£(k) [17] is as follows 


H = H ( o) + H (1) 


E(+) 0 

0 -E ( _) 


S R 
R+ Q ’ 


(A3) 


where E (+ ),E ( _), S, R and Q are certain N x N —dimensional submatrices. Systems 
underlying the matrix H and details of its construction (see e.g.[15,17]) are of no 
importance here. Let us note only that the relevant 2 N —dimensional basis set {T} is 
assumed to consist of two well-separated A - —dimensional subsets {'P( + )} and {T(_)}. 
The minus sign in front of E ( _) of Eq.(A3) is introduced for convenience. The super¬ 
script + designates the transposed (Hermitian- conjugate) matrix. 

The above-exhibited form of the initial Hamiltonian matrix H allowed us to look 
for the CBO matrix P, separate members of the power series of which (P(fc)) also 
are divisible into four submatrices (blocks). Moreover, a new version of the Raylcigh- 
Schrodinger perturbation theory (RSPT) has been formulated, wherein entire subma¬ 
trices (blocks) of the matrix H (i.e. non-commutative quantities) play the central role 
instead of usual (commutative) matrix elements. Accordingly, the new PT [19,20] has 
been called non-commutative RSPT (NCRSPT). As a result of its application, the 
corrections P^ take the following form 


p (fc) 


X(fc)+ G(fc) 

%) -x (fc) - 


(A4) 


where G^) are the so-called principal matrices of the NCRSPT determined by cer¬ 
tain matrix equations [18-20]. Meanwhile, diagonal positions of the corrections P^ 
are occupied by matrices X(*.) + and X(/q_ that have been referred to as intrasubset 
population matrices. These correspondingly refer to subsets {T( + )} and {H7(_)} and 
are expressible in the form of sums of products of matrices G(m of lower orders as ex¬ 
emplified below by Eqs. (A9) and (A15). The original derivation of members (P(k)) 
of power series for the matrix P [18] embraced terms to within second order only. 
Nevertheless, it is easily extendable to any k. 



As already mentioned, we confine ourselves here to a particular case of the matrix 
H of Eq.(A3). To this end, let us accept the equality E(+)= E(_) = I. [Besides, the 
resulting simplified Hamiltonian matrix coincides with that following from the initial 
matrix of our study shown in Eq.(l) after transforming the latter into the basis of 
bond orbitals (BOs) { 93 } [15]. (submatrices I and —I correspondingly represent one- 
electron energies of bonding BOs (BBOs) and of antibonding BOs (ABOs)]. The 
most important advantage of the above-specified condition consists in the possibility 
of an algebraic solution of matrix equations determining the principal matrices of the 
NCRSPT G(k), k = 1,2, 3.... As a result, these matrices meet the following recurrence 
relations 

G(fc) = -^|(SG( fc -i) - G( fc -i)Q) - L (fc) , (A5) 

where L(fc) are products of the same matrices of lower orders, e. g. 

L(i) =L( 2 ) = 0, L (3) = 2G(i)G^G(i), 

L( 4 ) =G( 2 )G^G(i) + 2 G(i)G^G(i) + G(i)G^G( 2 ) ,etc. (A 6 ) 

Finally, a useful relation R = —2Gp) follows for first order matrices in this case. 
Substituting the latter relation along with Eqs.(A3) and (A4) into Eq.(A 2 ) yields the 
following expressions for separate components of the sixth order energy, viz. 

4 “/ = - 47V(X (6)+ ), (A7) 

4? = - 2Tr(X, 5)+ S - X( 5 ,_Q)+ 8 rr(G( 5 )G+ ) ). (A 8 ) 

Let us now consider the components of Eqs.(A7) and (A 8 ) separately. The matrix 
X (6)+ determining the first component £^j is expressible as follows 

x (6)+ =G(5)G( i " 1) + G ( i)Gj } + G(4)GJ } + G(2)GJ } + G(3)GJ } + 

+ SG^G^G^G^G^G^) + G(i)G^ ) G( 1 )G^ ) + 

+ G(i)G^ ) G( 3 )G^) + G^jGJjG^jGJ) + G(3)Gj ) G(i)G[ l 1) + 

+ G( 1 )G^)G(2)G^ ) + G(2)G^ ) G(i)G^ ) + G(i)G^ ) G(i ) G( 2) + 

+ G( 2 )GJ ) G( 1 )GJ ) + G( 1 )GJ ) G( 2 )GJ ) + G( 2 )Gj ) G( 2 )G( l 1) . (A9) 

Substituting the above formula into Eq.(A7) shows that the component £^ contains 
matrices G(^) up to k — 5. The next step then consists in eliminating the matrices 



G( 5 ) and G( 4 ) from the expression concerned on the basis of our previous experi¬ 
ence when dealing with similar relations. To this end, let us take TV(G( 5 )G^) and 
Tr(G( 4 )GiL) separately and substitute recurrence relations of Eq.(A5) for G( 5 ) and 
G( 4 ), respectively. Moreover, Eq.(A6) also should be used along with allowed cyclic 
transpositions of matrices inside the Trace signs. The results of these procedures are 
as follows 

Tr(G (5 ) G+) =Tr(G {3) G+ 3) ) - 2Tr(G (3 ) G+G (1) G+) - 2Tr(G (1 ) G+G (2) G+) 

- 2rr(G+G (1) G+G (2) ) - ■2Tr(G (1 ) G+G (1) G+) 

- 2Tr(G (1) G+ 1) G (1 ) G+G (1) G+), (A10) 


Tr(G (4 )G+) =rr(G ( 3 )G+) + 2Tr(G ( 3 )G+G (1) G+) -Tr(G (1) G+G (2) G+) 

- Tr(G+G (1) G+G (2) ) - 2Tr(G (1 ) G+G (1) G+) (All) 

Employment of Eqs.(A7) and (A9)-(All) then yields the following formula for 

= - 20Tr(G (3 )G+) - 16Tr(G (3) G+G (1) G+) + 16Tr(G (1) G+G (2) G+)+ 

+ 16Tr(G+G (1) G+G (2) ) + 24Tr(G (1) G+G (1) G+) 

+ 8Tr(G (1) G+ 1) G (1) G+ ) G (1 ) G+) (A12) 

that, in turn, may be simplified considerably after eliminating Tr(G( 3 )G^G(i)G^). 
To this end, the following relation should be used, viz. 


Tr( G (3) GJ ) G (1 ) G+) =rr(G (1 ) G+G (2) G+) - 2Tr(G (1) G+ 1) G (1 ) G+G (1 ) Gj ) ) 

+ Tr(G+G (1) G+G (2) ) - Tr(G (1 ) G+G (1) G+). (A13) 

Derivation of Eq.(Al3) may be carried out analogously to those of Eqs.(AlO) and 
(All). The final expression for £^ is then as follows 

f ( < “ > = -20Tr(G (3 ,G+ ) )+40Tr(G ( i ) G+ ) G (1) G+ ) )+40Tr(G( 1) G+ ) G (1) G+ | G (1) G+ | ). 

(A 14) 

The second component of the sixth order energy (£^) also may be reformulated 
similarly. The first step of the relevant procedure consists in substituting into Eq.(A8) 
the expressions for X( 5 )+ and X( 5 )_, viz. 


X(5)+ - G (4)G("i) + G (i) G J) + GfsjGjj + G( 2 )GJ } + G(i ) Gj ) G ( i)Gj ) + 



G(i)Gj)G( 2 )Gj) + G(1 ) GJ ) G(1)G^ 1) + G^jGJjG^GJ), (A15) 
X (5)— = G 5)G(1) + G (1) G (4) + G (3) G (2) + G (2) G (3) + G (l) G (l) G (l) G (2) "f" 

G (l) G (!) G (2) G ( x ) + + G (1) G (2) G (1) G (1) + G (2) G (l) G (l) G (!)' ( Al6 ) 

Thereupon, relations of Eqs.(AlO), (All) and (A13) should be invoked to eliminate 
TV(G( 5 )GtA Tr(G( 4 )G^) and Tr(G( 3 )Gt,G(i)G^), respectively. The result is as 
follows 

£<« = 24^(0(3)05,)-48rr(G (1| G+ ) G (1) G+ ) )-48Tr(G (1) G+ ) G (1) G+ ) G (1) G+ ) ). 

(A17) 

After summing up the two components of the sixth order energy shown in Eqs. (A14) 
and (A17) in accordance with Eq.(A2), we finally obtain 

£ m = 4Tr(G, 3) G+,) - 8Tr(G (1) G+ ) G (1 ,G+ ) ) - 8rr(G (1) G+ | G (1) G+ 1) G (1) G+ 1) ). 

(A 18) 

The expression of Eq.(18) seems to be the most compact form of the correction 
concerned. However, it is not the most convenient one for practical applications, 
especially for the attempts of finding relations between separate increments of the 
overall correction £(@), on the one hand, and individual interorbital interactions, on 
the other hand. The main reason for that consists in the rather involved nature of the 
matrix G ( 3 ). Indeed, the relevant definition shown in Eqs.(A5) and (A6) embraces 
the product G(i)Gt,G(i) (along with the G( 2 )— containing term) that gives birth to 
increments like that of the last term of Eq.(Al8). To be able to sum up these similar 
increments, let us define a new matrix G° 3 ^ coinciding with the G( 2 )— containing term 
of Eq.(A5) for k = 3 as shown by the first relation of Eq.(12). Accordingly, the matrix 
G( 3 ) of Eq.(Al8) may be replaced by G^ — 2G(i)G^G(i). Thereupon, we may get 
rid of the newly-emerging term Tr(G° 3 ^G^G(i)G^) by constructing a relation like 
that of Eq.(Al3). The final formula for £^ is then as follows 


%) =4Tr(GfoGg) + 8rr(G (1) Gj ) G (1) G+ ) G (1) G+) - 16Tr(G (1) G+G (2) G+ 2) ) 

- 16Tr(G+G (1) G+G (2 )) + 8Tr(G (1) G+G (1) G+). (A19) 



Separate terms of the above expression are exhibited in Eqs. (8)—(11) and discussed 
nearby. It deserves adding here that third and fourth increments of Eq.(Al9) prove to 
be uniform in the case of AHs owing to the skew-symmetric nature of matrices Gp) 
and G( 2 ) [28]. After summing up these increments, a single destabilizing component 
of the sixth order energy (£^) arises (see Eq.(10)). 
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